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ABSTRACT 

Two classes of Gamma Ray Bursts have been identified 
so far, characterized by T 90 durations shorter and longer than 
approximately 2 seconds. We show here that the BATSE 3B 
data allow a good fit with three Gaussian distributions in logT 90 . 
The x 2 statistic indicates a 40 % probability for two Gaussians, 
whereas the three-Gaussian fit probability is 98 %. Using another 
statistical method, it is argued that the probability that the third 
class is a random fluctuation is less than 0.02 %. 

Subject headings: gamma rays: bursts 
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1. Introduction 

In the BATSE 3B catalog ( Mccgan, et al. 1996 ) there are 1122 Gamma- 



Ray Bursts (GRBs), of which 834 have duration information. [Kouveliotou 
|et al. 1993| have identified two types of GRB based on durations, for which 
the value of T 90 (the time during which 90% of the fluence is accumulated) 
is respectively smaller or larger than 2 s. This bimodal distribution has 
been further quantified in another paper ( |Kouveliotou, et al. 1995] ), where a 



two-Gaussian fit is made. In order to make further progress in quantifying this 
classification, one of the issues which needs to be addressed is an evaluation 
of the probabilities associated with a bimodal, or in general multimodal 
distribution. In this paper we take a first attempt at this, evaluating the 
probability that the two populations are independent, and we consider in 
addition whether a third group of bursts can be identified. 



2. GAUSSIAN FITS IN logT 90 

For this investigation we have used a smaller set of 797 burst durations 
in the 3B catalog, because these have peak flux information as well, and we 
use the Tg measures provided in this data set. For a one- Gaussian fit the 
X 2 probability is less than 0.1 %, so the one-Gaussian null hypothesis can 
be rejected. A log-normal two-Gaussian fit using 52 time bins with bin size 
0.1 in the log (fig. 1.) has a 40 % probability. This indicates that this null 
hypothesis cannot be rejected ( [Press, et al. 199"2]) . However, in the middle we 
have three consecutive bins, which have extremely large deviations. 

The interesting three bins in the middle have 3.2 a, 4 a and 2.5 a 
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deviations (fig. 2.). The probability to get such deviations in 52 bins is 6.3 % 
, 0.3 % and 58 % . Because we obtained these in three consecutive bins the 
probability that these deviations are random is ~ 10~ 4 . 

Such a deviation is interesting if it is in one/two/etc. bins. This raises 
the probability approximately by a factor five. These three bins contain 64 
GRBs. The two-Gaussian fit would account for 19, leaving 45 extra. The 
difference between these two numbers is more than 5 a, which corresponds to 
a probability of less than 6 x 10~ 7 . We can ask what is the probability of an 
excess in n (n small) consecutive bins located somewhere else. The number of 
trials is of order the number of possible locations of excess times the number 
of possible widths of excess, which is roughly 50 x 6 = 300. Therefore if the 
GRB duration distributions are log-normal then the probability that there is 
a deviation in these three bins from a two-Gaussian fit is more than 99.98 %. 

Because the duration distribution of the previously known two classes of 
GRB are nearly log-normal, we made a fit including a third Gaussian. The 
three-Gaussian fit (with nine parameters) has a x 2 = 24, which implies a 98 
% probability (fig. 3.). Table 1. contains the parameters of this fit. Looking 
at the goodness of the three- Gaussian fit, one might conclude that there is an 



"intermediate" type of GRB. This is an agreement with the |Mukherjee, et aL 



1998| result, who used a multivariate analysis and find that the probability of 



existence of two clusters rather than three is less than 10 4 . 

Both the two- and three-Gaussian fits have acceptable x 2 values, 
therefore neither can be rejected as an inappropriate fit. If we view the 
problem as one of model comparison, then, if the third Gaussian does not 
exist, the change in x 2 between the two fits from adding the three-parameter 
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third Gaussian should be distributed as x 2 °f three degrees of freedom ( |Band 



et al. 1997|) . The change in \ 2 °f 46.8-24.0=22.8 implies a chance probability 
of less than 10~ 4 indicating that a three Gaussian fit is a highly significant 
improvement over a two Gaussian fit. 

If we had considered six bins in the "intermediate" region, the two- 
Gaussian fit would show an excess in these bins of 61 GRBs. We may estimate 
that the "intermediate" group contains more than 50 but less than 70 GRBs, 
which would represent ~ 8 % of the GRB population. 



3. SYSTEMATICS 

The BATSE on-board software tests for the existence of bursts by 
comparing the count rates to the threshold levels for three separate time 
intervals: 64 ms, 256 ms, and 1024 ms. The efficiency changes in the region 
of the middle area because the 1024 ms trigger is becoming less sensitive 
as burst durations fall below about one second. This means that at the 
"intermediate" timescale a large systematic deviation is possible. To reduce 
the effects of trigger systematics in this region we truncated the dataset to 
include only GRBs that would have triggered BATSE on the 64 ms timescale. 

Using the Current BATSE catalog CmaxCmin table ( [Meegan, et al. 
|1998| ) we choose the GRBs, which numbers larger than one in the second 
column (64 ms scale maximum counts divided by the threshold count rate). 
Because this process reduced the bursts numbers very much we used the 
4B catalog ([Paciesas, et al. 1998|) data. In the 4B catalog there are 1234 
GRBs which have duration information; unfortunately just 605 burst satisfied 
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the above condition. In the interesting three bins we have 35 GRBs. The 
two-Gaussian fit would account for 13, leaving 22 extra. Therefore 63 % of 
the population are still a "deviation" (note in Sect. 2. it was 70%). 

This is still a 4-5 a deviation, depending on whether one uses the 
expected number or the observed one. Therefore after neglecting some 
systematics the probability that there is a deviation in these three bins from 
a two-Gaussian fit is about 99.8 % or 99.98 %. 



4. LOG N - LOG P 

In Fig. 4 we show the log iV — log P distribution of the two previous 
classes of GRB (see also |Kouveliotou, et al. 1995| and [Horvath, et al. 1996| ) . 



and compare this with the log N — log P in Fig. 5 for the "intermediate" 
class of bursts, consisting of 74 events (logT 90 =0.4-0.8). The logN-logP 
distributions of the T 90 <2.5 sec and the Tg >6.3 sec are on the fig. 4. 
The former has a Euclidean part, and the left portion goes over into a -1.05 
slope. The second group has a -1.06 slope. This corresponds to the results of 
[Kouveliotou, et al. 1995 . 



The logN-logP distribution of the "intermediate" type GRBs (fig. 5.) 
may be different, but this is questionable. 



5. HARDNESS RATIOS 



The average hardness ratios of the short and long bursts are 4.5 and 2.5. 
The "intermediate" type has 2.2 average hardness ratio. If this group were a 
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random mix of some shorter long-type bursts and some longer short-type of 
bursts, the expected value should have been between the two average hardness 
ratios. The actual value differs significantly from this, being close (and even 
somewhat smaller) than the hardness ratio of the long bursts. However, from 
a purely hardness-ratio point of view, it does not seem possible to distinguish 
the "intermediate" bursts from the long ones, thus the hardness information 
does not support, but also does not contradict, the "intermediate" burst class 
hypothesis. 



6. CONCLUSION 

It is possible that the three log-normal fit is accidental, and that there 
are only two types of GRB. However, if the Tg distribution of these two types 
of GRBs is log-normal, then the probability that the third group of GRBs is 
an accidental fluctuation is less than 0.02 %. The logN-logP distribution and 
hardness ratio are also suggestive of the fact that the intermediate duration 
(T 90 =2.5-7 sec) bursts represent a third class of GRBs. 

This research was supported in part through NASA5-2857, OTKA 
T14304, Szechenyi Fellowship and KOSEF. Useful discussions with E. 
Feigelson, E. Fenimore, A. Meszaros, P. Meszaros, and J. Nousek are 
appreciated. The author also thanks an anonymous referee for useful 
comments that improved the paper. 



- 8- 



REFERENCES 

Band, D. L, et al. 1997, ApJ, 485, 747 

Horvath, L, Meszaros, P., & Meszaros, A. 1996, ApJ, 470, 56 
Kouveliotou, C, et al. 1993, ApJ, 413, L101 

Kouveliotou, C, et al. 1995, Third Huntsville Symposium on GRB. New 
York: in AIP Conference Proceedings 384, 84-89. 

Meegan C. A., et al. 1996, ApJS, 106, 65 

Meegan C. A., et al. 1998, Current BATSE Gamma-Ray Burst Catalog, on 
the Internet [ht t p : / / www . bat se . msf c . nasa . gov / dat a/ gr b / cat alog/ 



Mukherjee, S., Feigelson, E. D., Babu, G. J., Murtagh, F., Fraley, C, & 
Raftery, A. Submitted to ApJ. |astro-ph/9802085 



Paciesas, W. S., et al. 1998, The Fourth BATSE Gamma- Ray Burst Catalog. 
ApJS, in preparation. 

Press, W. H., Teukolsky, S. A., Vetterling, W.T., & Flannery B. P. 1992, 
Numerical Recipes in Fortran, Second Edition, Cambridge University 
Press, Cambridge 



This manuscript was prepared with the AAS IATjrjX macros v4.0. 



- 9- 



Figure caption 1. 

Distribution of log T 90 for 797 bursts from the 3B catalog. The solid line 
represents a fit of two log-normal Gaussians using 6 parameters and 52 bins. 
The best fit x 2= 46.8, which implies a 40 % probability. 

Figure caption 2. 

Difference between the two log-normal Gaussian fit and the logT 90 
distribution. The solid lines represent one sigma errors of the T 90 bins. 
Dotted lines and dashed lines mean two and three sigmas. Except for the 
three middle bins, only one bin has ^ 2 sigma deviation, and seven bins have 
i<5 1 sigma. 

Figure caption 3. 

log T 90 distribution. The solid line represents a fit of three log-normals. 
The x 2 value implies a probability of 98 % . 

Figure caption 4. 

The logN-logP distributions of the two type of GRBs. The left curve is 
for T 90 <2.5 sec (223 GRBs), while the right curve is for T 90 >6.3 sec (508 
GRBs). 

Figure caption 5. 

logN-logP distribution of the "intermediate" type (2.5<T 90 <6.3) GRBs. 
The slope of -0.65 is substantially different from that of the two other types 
of GRBs. This is compatible with the "intermediate" GRB being a separate 
class of objects. 
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Table 1. 





raiding s) 


a(lgs) 


memb 


% 


^90 


first 


- 0.35 


0.50 


236 


30 


< 5 s 


second 


1.52 


0.37 


497 


62 


> 4 s 


third 


0.64 


0.14 


61 


8 


2.3-8 



The parameters of the three- Gaussian fit of GRBs. 

Average duration (mid), width of the Gaussian (a) (both in lg sec). 
Members of the group (memb) and population of the group (%). Border of 
the group in T 90 (unfortunately they overlap). 



